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SUMMARY 

A simple interaction curve is presented for evaluat- 
ing the conditions of combined shear and direct stress 
under which an infinitely long, flat plate with equal 
elaetio restraints against rotation along the edges will 
become unstable. The theoretical work that led to the 
interaction curve is presented in the form of appendixes. 

INTRODUCTION . 

In the design of stressed-skin structures, considera- 
tion must sometimes be given to. the critical etresees for 
a sheet under a combination of shear and direct stress. 
The upper surface of a wing in flighty. for example, may 
be subjected to combined shear and compressive stress and 
the lower surface to combined shear and tensile stress 
simultaneously. Under this condition the upper surface 
may buckle at a lower compressive stress than if the shear 
were not present and the critical shear stress for the 
lower surface will be increased by the presenoe of the 
tensile street?. The purpose of the present report is to 
specify the combination of shear and direct stress under 
which such, surfaces may be expected to buckle. 

INTERACTION CURVE 

f 

For the case of flat sheet between identical stlff- 
enere, the combination of shear etress T and compres- 
sive stress - o$ that will cause buckles to appear in the 

sheet has been exactly solved in appendix A. The values 
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of ff x and T may oe given In the form (reference 1, 
pp, 366 and 369) 



where 
D 



u 

b 
t 

k. 



T 3 



b a t 



a 

k o TT D 
b a t 



flexural- stiffness of sheet per unit length 
St 3 



13(1-H 3 ) 

modulus of elasticity 
Foi«8on 1 s ratio 

wiflLth of sheet "between stiffeners 

thiokness of sheet 

( ooef f ici Put a dependent upon the restraint supplied 
* "by the stiffeners 

Appendix A describes the intricate relationship that 



exists between k. 



and 



As is usual with exact so- 



lutions of buckling problems, this solution is very labo- 
rious to apply. 



The possibility of great simplification in the com- 
putation of the critical stresses is indioated by the re- 
sults of the application of an energy method presented in 
appendix B. In this appendix it is shown that if the re- 
straint supplied to the sheet by the stiffeners is inde- 
pendent of the wave length of the buckles, the following 
relation exists: 

R q + E 8 3 = 1 (1) 



where 
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H # ratio of shear stress when buckling occurs In combined 
shear and direct stress to shear stress when "buck- 
ling ooours In pure shear 

H ratio of compressive stress when buckling oocurs In 
combined shear and direct stress to compressive 
etreee when buckling oocurs In pure compression 

Inasmuch as design charts for determining the critical 
etresses under pure shear and pure compression for flat 
sheet with equal elastic edge restraints have been pre- 
sented in references 2 and 3, equation (1) may be easily 
applied to the combined loading of such sheet. In an ac- 
tual structure, however, the restraint Is usually depend- 
ent upon the wave length of the buckles, and the formula 
must therefore be shown to hold with sufficient accuracy 
In this case also. 

Bquation (1) was tested as described in appendix 0, 
by use of exact values of the critical stresses obtained 
by the method of appendix A, The results are shown In 
figure 1, in which the curve represents the values ob- 
tained by the combined loading formula of equation (1) 
and the plotted points are for the exact values oomputed 
by the method of appendix A and listed in table 1. In 
some cases, the restraint was assumed to be independent 
of the wave length; whereas, in others, a relationship be- 
tween restraint and wave length typical of a sturdy stiff- 
ener (referenoe 4) was assumed. Some values were also 
computed for tension in the sheet (k 0 negative). .All 

values are seen to lie on the curve or close. to it, indi- 
cating that equation (1) holds with sufficient accuracy 
for engineering purposes in all eases. 

00N0LUSI0JTS 



It is concluded that the values of combined shear and 
direct stress at which an Infinitely long, flat plate sup- 
ported along the edges with equal elastic restraints 
against •rotation will become unstable may be determined 
for practical engineering purposes by the equation 

a 

H + H =1 
o s 
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where ■ 

E c ratio of direot stress when "buckling occurs in coin- 
Dined shear and direct stress to compressive stress 
when Duckling occurs in pure compression. Tension 
is regarded as negative compression. 

H B ratio of shear stress when buckling occurs in com- 
bined shear and direct stress- to shear stress when 
buokling occurs in pure shear 
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APPENDIX A 



SOLUTION BY DIJPFEHENTI AL EQUATION 



■The exact solution for the critical stress at which 
buckling ocours in a flat rectangular plate subjected' to 
combined shear and compression forces in its own plane 
may' be obtained'by solving the differential equation that 
expresses the equilibrium of the buckled plate. The plate 
is as.sumed to be infinitely long, and equal elasti.o. re- 
straints against rotation are assumed to be present along 
the two edges of the plate. 

Tigure 2 shows the coordinate system used. The dif- 
ferential equation for equilibrium of a flat plate under 
shear and direct stress in the direotion of one axis is 
(reference 1, p. 306) 



4 

d w 



+ a 



0 w 



ox* dx 8 dy a 



+ 3Tt 



.3 

0 w 
dxdy 



+ ff T t 



bx a 



= 0 



(Al) 



where 

D flexural stiffness of plate per unit length 
3 modulus of elasticity 



Bt 3 1 
12(l-u a )J 
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ro 
i 

1-3 



H Poiaaon^ ratio 
t thickness of . plate 

w deflection of plate at (x,y) from unstressed position 
T uniformly distributed Bhear stress 

cr x uniformly distributed compreeelTe stress in the direc- 
tion of x 

It ie known that the oritioal stresses are of the form 
(referenoe 1, pp. 356 and 359) 



k 0 '-" a D 
b a t 



lea" D 
b a t 



(A3) 



are 



where b is the width of the plate and k c and k 

constants which, for an infinitely long plate, depend only 
upon the elastic restraints along the edges of the plate. 
Substitution of expfepsions (A2) in equation (Al) gives 



^4 st ~4 _ a, v s a, ~b 

^ + 3-V^ + M + £!L T £iL — -^H" 0 (A3) 
dx* dx s oy 8 by 4 b B &xdy Lb dx 



If the plate la infinitely long in the x-direotion, 
all displacements must be periodic in x, and the defleo- 
tion surface may be taken in the form 

w = Ye * (A4) 



where 7 is a function of y only and ^ .is the half 
wave length' of the buokl'efs in the x-direotion. 

Substitution of expression (A4) into the differential 
equation (A3) gives as the ecuation which determines 7 
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d T 2n d T it 2tt lk B dY n k fl 

— + — T + - T «= 0 (A6) 

dy \ a dy 8 X 4 b°* dy b S X a 



A solution of equation (i.5) ia 

T = e 

when m ia a root of the characteri atio equation 

4 ^/^V* 8 a/Vb\ /VbV aAr-bV „ . v 

m + 2 w ■ - ^ W* 8 " + W " W w k ° = 0 (A6> 

Except for the aubatitution of - tt 3 (^ k c f or 

equation (A6) ia identical with equation (A6) of reference 
2, in which equation (Al) of this appendix was solved with 
k Q = 0. Vith thie change, all the results obtained in 

that appendix are applicable here. The stability crite- 
rion for combined compression and shear is therefore for- 
mally the same as for shear alone as given by equation 
(A19) of reference 2. This criterion is 

2aP^4Y 3 - ~ ^^cosh 2a cos 2P - cos 4Y^ 

r 3/ a a\ / a aV / a a a\ € a l 
- |4Y - a J -(^p + a j ~ (4Y - p + a J± Jeinh 2a Bin 20 

+ € J^a^Y 3 + of+ 0 3 ^ cosh 2a sin 2g 

+ p^4Y 3 - a 3 - P 3 ^ Binh 2a cos 2P - 4aPY 8 in 4Y^ = 0 (A7) 



where the relation between k B and a, p, and Y is 



the same ae in reference 3, equation (A23) 



"but In the present report 



r (A9) 



The restraint coefficient e in a measure of the relative 
resistance to rotation of the restraining element at the 
edge of the plate; it is discussed more fully in refer- 
ence 3. 

The procedure for evaluating k 8 , after values of 
k 0 and an expression for £ have been chosen, is as 

follovs: A value of "b/X is assumed; a series of values 
of y are taken until one is found which, together with 
the corresponding values of a and 3 as computed from 
equations (A9), satisfies equation (A7); k B is then com- 
puted from equation (A8). Another value of b/A is as- 
sumed, a new Bet of values of Y, a, and 3 are found 
that satisfy the stability oriterion, and a new value of 
k a is computed. The entire process is repeated until by 

plotting k a against b/X the minimum value of k 8 is 

found. In the case where e is a function of b/X, e 
must be re-evaluated each time a different value of b/X 
is assumed. This minimum value of k B and the chosen 

value of k 0 , when inserted in equations (A3), give a 

critical combination of shear and direct stress. 



8 



APPENDIX B 



SOLUTION BY ENERGY METHOD FOR CONSTANT t 



In appendix B of referenoe 3, an energy solution was 
given for the stability of a plate under pure shear by 
use of oblique coordinates . This solution may be extended 
to cover the oase of combined shear and' compression load- 
ing by simply adding to the work done by the shear force 
an additional term expressing the work done by the compres- 
sive force. 

Figure 3 shows the coordinate system and the plate 
dimensions. The procedure is to evaluate the terms in the 
e qua ti on 

T c + T e ■ V i + V a (B1 > 

where 

T Q work done by the compressive force per half wave length 

T 8 work done by the shear force per half wave length 

Vj strain energy in the plate per half wave length 

T a strain energy per half wave length in the elastic 

restraining members assumed to be present along the 
edges of the infinitely long plate 

The deflection surface is taken to be the same as in 
reference 2, equation (B2); that is, 



w = w 0 



[¥(£ - § ♦ ( l + ?>» 5f] «■ V <»> 



where € is the restraint coefficient defined as in ref- 
erence 2. The valued of T 8 , 7 lt and V a may be taken 

directly from eauations (B3), (B4), and (B5) of reference 
2: 

tt bjt sin <t> 



T s = w o 



2X 



IA120 8 ti / \2 Tt s 2j 



° 4D a x cob <j> LViao 8 W Va tt»/ sj 

*©" [g-*>'*g-*vs 

♦ 8(1 + a . ln ,*,[(i_.L) € . + (i_ JL)« + |]} 



T 3 . ... , (35) 

The work T c done "by the compressive stress o T is given 
hy 



By use of the assumed expression (B3) for v, T Q is 
found to he 

- - .■ ^(^* i- (|- C.) 

When equations (B3). (B4) , (B5) , and .(B6) are sub- 
stituted in eauation (31), the result is 

a x = k 0 ^ - 2t tan <D (B7) 



where 
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©■ . °> 



„ 1 + 2 



C08*<J> /V\ 3 C08 3 <|> 



Vl30 8 it 8 / \3 ttV 3 

(JL _ iV. + (1 _ 0) f + i 
\24 ttV \.3 it'/ 3 



V130 S iiy \^3 ttV 3 



The angle of Inclination <t> and the ratio b/X will 
adjust themaelvea to make o y a minimum; that is, 



* CT * - 0 (38) 



ocr x 



(B9) 



If the shear stress T is considered br s given constant, 
and the restraint coefficient e is independent of b/A, 
the values of b/A and <j> resulting from the operations 
indicated by eauations (38) and (39) are: 

' coa * (B10) 
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<{> tan" 1 



X 



(Bll) 



Inasmuch as the value of b/\ given by equation (BIO) 
and the value of <|>- given "by equation (Bll) are the val- 
ues vhich make the stress a x a minimum, it is necessary 

to substitute them in equation (B?) to obtain the critical 
compressive stress. When this substitution is made and 
the result expressed in nondimensional form, it is found 



In order to Identify the denominators in equation 
(Bis), it is necessary to refer to the en-ergy solutions 
for pure compression and for pure shear. "Trom equation 
(B14) of reference 3, it is f.ound that 



that 



(B13) 




where kp Q is the value of k Q for pure compression; or, 
e being assumed constant and kp 0 minimized with respect 
to b/*, ' . " . . 




(BIB) 



Trom equation (B6) in reference 3, it is found that 



sin ie$ 



1 




C 1 cos a <b 



+ 0 a (l + 3 sin»4>) 
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where kp B is the value of k fl for pure shear. If the 
value of e is assumed constant and kp B is minimized 
with respect to b/x and <fr , .kp 8 becomes 



or 



k ps 8 = 



kp B - 2 + 2 /~C7 0 a + | 0 a a 

4^ + 3/07 0 a + | 0 3 a ^ = ^2yo7+ 3C^2/C7+ (B14) 



Substitution of relations (Bis) and (B14) in eauation 
(B12) gives the combined loading formula 



-i_ ( t Y 

k k "* 

*po K ps 



Or 



r~ + (ir) = 1 (B15> 

K pc ^ K ps' 

where k Q and k„ are written for and 

TT*D/b 3 t 

, respectively. Equation (B15) may also be wrlt- 

TT 3 D/b 3 t 

ten in the form 

a o + B s S = 1 ( Bee equation (l)) 

where 



a 0 



_k 
k 



8 k 



pc 

*S 



ps 
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These definitions of - Eg and - R a are. equivalent to 

the definitions given in terms of stresses under sections 
Interaction Curve and Conolusions in the "body of this re- 



The theory of appendix B, although "baaed upon an ap- 
proximate deflection surface for the Buckled plate, indi- 
cates a simple algebraic relationship between k c and 

k e , when the restraint coefficient e is independent of 

the wave length X, of the form 



k^g value of k Q when k B «= 0 (pure compression) 
kp 8 value of k a when k Q = 0 (pure shear) 

The exact solution found in appendix A provides the 
means by which equation (Cl) may he tested. An exaot nu- 
merical check could not be expected because some error 
would he introduced in reading the charte for kp 0 and 

kp B and in the graphical work involved in obtaining k fl ; 

the test should reveal", however, whether' the formula is 
sufficiently accurate for engineering computations. 

Equation (Cl) was tested for the following values of 

e : 



port . 



APPENDIX C 



TEST Off INTERACTION 70EKULA 




(Cl) 



where 



c 



10 




14 



The first three values of e are independent of the 
wave length, and are therefore representative of the val- 
ues of £ upon whioh equation (CI) ie baaed. The final 
value of e varies with the wave length in a manner typ- 
ical of a sturdy stiff ener (reference 4), and this varia- 
tion is therefore representative of a practioal oase. 

Tor eaoh value of € a value of k Q was chosen, 

and the associated value of k Q required to cause the ' 

plate to buokle was computed as described in appendix A, 
lour oases were computed with tension on the plate (k Q 

negative). The values of k pc and kp fl were then read 

from their respective charts in references 2 and 3. Tor 
the oase of oonstant e, the values of k pc and kp B 

were read at the minimums of their appropriate e-ourves. 
(The values of b /\ at which the plate buckles under pure 
compression, pure shear, and under the oombined loading 
are, of course, all different.) Tor' the oase of variable 
€', it is necessary to shift from one e-ourve to another, 
in conformity with the assumed relationship between e 
and X, until minimum values of kp. c and kp B are found. 

The results of the numerical test of the interaction 
formula are shown in table I and in figure 1. . The final 
column of table I indicates the value of the left-hand mem- 
bers of equation (Cl). Tor the cases in which e 1b in- 
dependent of X the difference in this value from unity 
is only a fraction of a percent, except when H c has a 
large negative value, in which case the difference 1b 8 to 
5 peroent; the values in the last column of table 1 will be 
more in error for these values of R 0 because their com- . 
putation involves subtracting quantities of the same order 
of magnitude. Tor the oase in which e varies with X 
the difference from unity is 1 or 2 peroent. The validity 
of the interaction formula is therefore established for all 
engineering purposes. 
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TAB LIS I 

TBI ST 07 ACCURACY 0! INTERACTION FORMULA 
[Values oomputed hy the method of appendix A] 



k 0 


k pc 


v. k ° 


*e 








£ b 0 


4.00 
3,00 
2.10 
1.00 
0 

-30.00 


4.00 


1.0000 • 
.7500 
.5250 
.2500 

0 

-6.0000 


0 

2.682 
3.690 
4.623 
6.335 
13.997 


6.335 


0 

.5027 
.6917 
.8665 
1.0000 
3.4362 


1.0000 
1.0027 
1.0034 
l.OOOB . 
1.0000 
.9468 


€ = 10 


6.605 

5.00 

2.80 

1.30 

0 

-1.30 
-6.00 


5.605 


1.0000 
.8931 
.4996 
.2319 

0 

-.2319 
-1.J0705 


0 

2.380 
6.115 
6.357 
7. 270 
8.052 
10.418 


7.270 


0 

.3274 
.7073 
.8744 
1.0000 
1.1076 
1.4330 


1.0000 
.9993 
.9999 
.9966 

1.0000 
.9949 
.9830 




6.98 
5.40 
3.49 
1.80 
0 

-18.00 


6.98 


1.0000 
.7736 
.5000 
.2579 

0 

-2.5788 


0 

4.272 
6.347 
7.737 
8.980 
15.953 


8.98 


0 

.4757 
.7068 
.8616 
1.0000 
1.8879 


1.0000 
.9999 
/9996 
1>. 0002 
I. 0000 
. 9854 
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6.119 
' 5.60 
3.00 
1.00 
0 


6.119 


1.0000 
.8998 
.4903 
.1634 

0 


0 

2.315 
5.018 
6.328 
6.883 


6.883 


0 

.3363 
.7290 
. 9194 
1.0000 


l.oq'oo 

1.0119 
1.0317 
1.0087 
1.0000 



1.2 
1.2 
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(a) Shear and 
compressi-h 



(o) Shear and 
tension 



Figure 1,- 

Interaction 

curve for 

plates 

under 

combined 

shear and 

direct 

stress. 
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Tigure 2,- Inf init&ly long rectangular plate tinder combined 
shear and compression; coordinate system used in 

appendix A. 
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Figure 3.- Oblique coordinate system used in appendix B. 



3 1176 01345 5424 



